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Abstract 

We show that solutions u{x,t) of the non-stationnary incompressible Navier-Stokes 
system in M.'^ {d > 2) starting from mild decaying data a behave as \x\ ^ oo as a potential 
field: 

where is a constant and Kh.k = /q (it/i|wfe)L2 is the energy matrix of the flow. 
We deduce that, for well locaHzed data, and for small t and large enough |a;|, 

ct|a;|-('^+i) < \u{x,t)\ < c't\x\-'^'^+^\ (ii) 

where the lower bound holds on the complementary of a set of directions, of arbitrary 
small measure on S''"^. We also obtain new lower bounds for the large time decay of the 
weighted-L^ norms, extending previous results of Schonbek, Miyakawa, Bae and Jin. 



Nouveaux profils asymptotiques 
des solutions non-stationnaires de Navier-Stokes 

On montre que la solution u{x,t) de I'equation de Navier-Stokes incompressible 
dans {d > 2) issue d'une donnee de Cauchy generique et moderement decroissante a 
se comporte, pour \x\ oo, comme un ecoulement potentiel donne par la formule (i) ; 
7d est une constante et KhM = Jo{uh\uk)L^ est la matrice d'energie de I'ecoulement. 

On en deduit que, si la donnee est bien localisee, le champ de vitesse verifie (ii) pour t 
suffisament petit et \x\ assez grand. La borne inferieure est valable sur le complementaire 
d'un ensemble de directions, de mesure arbitrairement petite dans S'^"^. On obtient aussi 
de nouvelles bornes inferieures du taux de decroissance en temps grand des moments de 
la solution dans qui etendent des resultats anterieurs de Schonbek, Miyakawa, Bae 
and Jin. 



Keywords: Asymptotic behavior at infinity. Decay. Upper bound estimates. Lower bound 
estimates. Mild solutions to the Navier-Stokes system. Asymptotic separation of variables. 
Peetre weight. 
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1 Introduction 



Let a be a divergence- free vector field in {d>2). We consider the Cauchy problem for the 
Navier-Stokes equations : 

dfU — An + (n • V)^ = — Vp, 
< div« = 0, (NS) 
^u{x, 0) = a{x). 

The unknowns are the velocity field u = {ui, . . . ,Ud) and the pressure p. The problem has to 
be solved on x [0,+oo) or at least on x [0,T) for some T > 0. 

Because of their parabolic nature, tlie Navier-Stokes equations feature an infinite-speed 
propagation effect in the space variable. This phenomenon is usually described by the fact that 
compactly supported initial data give rise to solutions which immediately have non-compact 
support. On the other hand, because of the pressure, which can be eliminated from the 
equations only by applying a non-local operator, the solutions of the Navier-Stokes equations 
have quite a different behavior as \x\ — ^ oo from that of solutions of non-linear heat equations. 

The main purpose of this paper is to study such asymptotic behavior. For example, we 
address the following problem. Assume that, at the beginning of the evolution, the fluid 
is at rest outside a bounded region (say, a G C^{W^), the space of smooth, solenoidal and 
compactly supported vector fields). At which velocity will the fluid particles that are situated 
far from that region start to move ? 

We will obtain sharp answers to this and related questions by constructing new asymptotic 
profiles of solutions to (NS), predicting the pointwise behavior of u as |x| — > -|-oo. 

A few asymptotic profiles of solutions to the Navier-Stokes equations in the whole space 
are, in fact, already known. For example, F. Planchon [22], studied self-similar profiles. 
However, his results cannot be used in the case of initial data decaying at infinity faster than 
\x\~^, since the only possible self-similar profile, in this case, would be the zero function. For 
faster decaying data, the asymptotic profiles of A. Carpio [6], Y. Fujigaki, T. Miyakawa [8], 
Miyakawa, Schonbek [21], T. Gallay, E. Wayne [11] and Cannone, He, Karch [5] provide 
valuable information about the large-time behavior of the velocity field. However, in all these 
works the asymptotics are obtained by computing some spatial norms of expressions involving 
the solution. The limitation of this approach is that most of the information on the pointwise 
behavior of the velocity field is lost. 

Our method is different, and consists in proving that, asymptotically, the flow behaves as 
a linear combination of functions of separate variables x and t. 

Our proflles imply that, without external forces, the flow associated with decaying initial 
data behaves at infinity as a potential field, with a generalized Bernouilli formula relating 
the pressure to the energy matrix {uh\uk)L^ of the flow. This illustrates the fact that the 
spatial behavior at inflnity of the flow is almost time-independent, contrary to the temporal 
asymptotic, which is known to be influenced by spatial decay. 

Notations 

1. We denote by the space of all measurable functions (or vector flelds) / on M*^, such 
that : 

||/|U^=ess sup(l + H)''|/(x)|<+oo. 
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The space ([0, T); L^) is made of functions u{x, t) such that u{t) G for all t > 
and 

' lim \\u{t') - u{t)\\L^ =0 if t > 0, 
u(t) u(0) in the distributional sense. 

2. For positive (f), the notation f{x,t) = Ot (^(a;)^-^) means that \4>{x)f{x,t)\ < Ct, for 
some function t ^ Ct, possibly growing as i ^ oo, but locally bounded. 

3. The solution of the heat equation is 

e*^a{x) = (47rt)-'^/2 / a{y) dy. 

4. We also adopt the standard Kronecker symbol: Sij = 1 ii i = j, and Sij = otherwise. 

Our starting point is the following well known result (see [19], chapter 25). 

Theorem 1.1 Let d > 2. There exists a constant 7 > such that for any divergence-free 
vector field a G L°°(M°'), one can find 

T > 7min {l; ||a||^^} 

and a unique mild solution u £ Cw ([0, T]; L°°) of (NS). This solution u is smooth for t > 0. 
Moreover, if a belongs to for some > 0, then we also have : 

ueCy, {[0,T];Lf) 

with ■& = minli?; d+1}. 

This conclusion can be restated in a slighlty different way (see also [24, Proposition 3]) : 

u{x, t) = e*^a + Ot ((1 + "1^2^;^+!}^ 

on [0,T] X M<^. 

Asymptotic behavior of local solutions 

We can now state our first main result. Let us introduce the energy matrices : 

rt 



£h,k{t) = (uhUk)iy,t)dy and Kh,kit) = I / {uhUk){y,s) dy ds. 

JRd Jo JRd 

The following theorem describes the asymptotic profile of local solutions. 



(1) 



Theorem 1.2 For '& > ^ and an initial datum a G Lf, let u G C^([0, T]; L^^) he the 
solution of (NS) given by the preceding theorem. The following profile holds for \x\ — > +00 ; 

u{x,t) = e*^a{x) + Vn{x,t)+Ot (^\x\- '^'''i^^ ' '^+^}^ (2) 

where Il{x,t) is given by : 

n(-'*) = 7^E(^-^)-^MW (3) 
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and 7d = TT '^/^ r(^i^). If, moreover, the first and second order derivatives of a belong to L"^, 
then there exists a constant po such that the following profile holds for t > : 

p{x,t)=po - 7.E - ^) ■ + 0. (|.|— (4) 

h,k 

Remark 1.3 This theorem essentially says that, for mild decaying data (this is the meaning 
of the assmnption a G L^, with -i? > ^^), 

u{x,t) ^ e*^a(x) + Vn(x,i), as \x\ — > oo. 

In other words, the flow behaves at infinity as the solution of the heat equation phis a potential 
field at infinity. In particular, if follows that for fast decaying data {i.e. when -& > d + 1, we 
simply have 

u{x,t) ~ Vn(x,t), as \x\ oo, 

since the linear evolution can be included in the lower order terms. 

Theorem 1.2 does not cover the case of slowly decaying data {i.e., the case ^? < ^). The 
spatial asymptotic of those slowly decaying solutions (inchiding self-similar solutions) has a 
different structure and cannot be constructed with the same method. We should consider it 
in an independent paper. 

Remark 1.4 The decay of the remainder in (2) cannot exceed jxl""^"^. Indeed, (NS) being 
invariant by translation, the choice of the origin is arbitrary and one can easily check that 

Vn(x-xo,i) - vn(x,t) 

decays at infinity as if IT ^ and xq 0. 

Even if u{x, t) develops a singularity in finite time, the potential field in (2) will remain 
uniformly bounded away from the origin : 

|Vn(x,t)| < C||a||^2t|a;|-'^-^ 

However, the above result provides no information about the singularity itself, nor does it 
prevent it from appearing : as long as the solution is smooth, the remainder of (2) compensates 
for the singularity at the origin of Vn(x,t). 

Remark 1.5 The above profile for the pressure has some analogy with Bernoulli's formula 
for potential flows : 

P = Po+ 2^^^- 

Such a formula holds rigorously for the stationary Euler equation with no external force, but 
this identity can also be useful when dealing with high Reynolds flows around aerodynamical 
bodies (see, e.g., the description of the Prandtl laminarity theory in [14, Chapter 9]). 

The asymptotic profiles of Theorem 1.2 are meaningful when the leading term does not 
vanish identically. It turns out that this is the case for generic solutions. Indeed, the next 
result provides a necessary and sufHcient condition for VII to be identically zero. 

Proposition 1.6 Let u as in Theorem 1.2 and K = {Kh,k)- For any t G [0, T], the homo- 
geneous function x i— > Vn(x, t) vanishes identically on if and only if the matrix K{t) is 
proportional to the identity matrix, i. e. 

yh,ke{l,...,d}, Kh,k{t)=a{t)6h,k (5) 

with a = \TTK. 
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This shows that VII does not vanish for generic flows. Conditions (5) also occur in the 
paper of T. Miyakawa and M. Schonbek [21]. It is shown therein that a high decay rate of the 
energy of the flow for large time is essentially equivalent to (5) holding in the limit t +oo. 

Such orthogonality relations can also be described in terms of vanishing moments of the 
vorticity lo = curl u of the flow (see, e.g. [11], [12]). Focusing on the vorticity, in fact, has 
crucial advantages in the study of the large time behavior of solutions, especially in the two- 
timensional case. We refer e.g. to the recent work of Gallay and Wayne [13] on the global 
stability results of vortex solutions. 

On the other hand, the large space behavior of the vorticity is less interesting than that 
of the velocity field. This can be shown by taking the curl(-) operator term-by-term in 
formula (2): the term curl(Vn) identically vanishes. The physical interpretation of this 
remark is the following: if we start with an initial datum with localized vorticity, then the 
vorticity will remain localized as far as the solution exists (this fact, of course, was already 
known) . 

In principle, it would be possible to extend formula (2) and to write a higher-order asymp- 
totic for u as \x\ oo. The above observation allows us to predict that all the higher-order 
terms of the expansion of the velocity field must be curl-free. Otherwise, there would be a 
limitation on the decay rate of lo as \x\ — > oo and this would contredict e.g. the results of [12], 
[18]. In other words, all the higher-order terms of the expansion oi u should be gradients. 



Large time asymptotics 

Under a suitable smallness assumption such as 

ess sup |x| |a(x)| < eq, 



one can take T = +oo in Theorem 1.1 (see, e.g., [20]). Moreover, the localization property of 
the flow persists uniformly. One has : 

|u(x,0| <C(l + |x|)-"(l+i)-''/2 

for any a,/3 > such that a + fj < minj-t? ; d + 1}. When i} = d ov d + 1, one needs the 
additional assumption that the above estimate already holds for e*^a (see [1]). In particular, 
these estimates imply that, for large t : 

||2 iC, ifi?>^ 

II^IIl^([o,],l^)<|^^^_,+^+.^ if^<^, 

for all e > (this bound also holds for e = 0, i!) ^ ^i^, but we will not use this fact). 
We can now give our asymptotic profile for global solutions. 

Theorem 1.7 Given ?? > let u{x,t) be a solution of (NS) on M+ x such that 

\u{x, t)\ < Co (1 + |x|)-" (1 + t)-^/2 

for any a,(3 >0 such that a + f3 < "d. Then, 

u{x,t) = e^^a{x) + VIl{x,t) + \x\-'^-^ e(^-^^= ; +n{x,t) (7) 

with the following estimates : 

|^(x,i)|<Ce-'=N'hlli2([o,,],i2) (8) 
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and, for any < a < min{l, i-}, and all t > 1, 

\n{x,t)\ < Ca|a:r'^-^-"r5+f , if^> d±3^ 

\n{x, t)\ < Ca,e kr<^-l-" t^-'^H-^ if ^ < ^ < ^. 

Due to the form of the remainder terms, it seems impossible to obtain a description of the 
pointwise behavior of u for large but fixed and t ^ oo. Conclusion (7) is interesting only 
for {x,t) such that |x| > C\/t + 1. For those points, this profile provides more information 
than those in [8] or [5] (on the other hand, our assumptions are necessarily more stringent). 



Applications 

For smooth and fast decaying initial data, according to Theorem 1.2 one has 

\u{x,t)\<Ct{l + \x\)-(''+'l (10) 

Theorem 1.2 allows us to answer the more subtle problem of the validity of the corresponding 
lower bound to (10). 

A first difficulty is the following: the upper bound ensures that u{-,t) is integrable, so that 
the divergence-free condition implies 

Vt>0, Vj G {l,...,d}, / Uj{x,t)dx = 0. 

jRd. 

In particular, since u is smooth for t > 0, no uniform lower bound by a given positive function 
can hold. However, non-uniform and anisotropic lower bounds do hold, even if the initial data 
is rapidly decreasing. 

More precisely, for generic flows (i.e. if we exclude flows with special symmetries) starting 
from fast decaying data, we will prove that for some small to > (depending only on the 
initial datum), and for j = 1, . . . , d, 

ct\x\-^'^+^^ < \uj{x,t)\ < c't\x\-^'^+^\ c,c' > (11) 

for all t G (0, io] and all > C/\/i, with x outside a small set of exceptional directions, 
along which the decay can be faster. In other words, the constant c in (11) is independent 
of t or |a;|, but does depend on the direction x/\x\ (see Theorem 3.1 below for a more precise 
statement). For example, we will see that in dimension two the exceptional set is made of 
at most six directions. The remarkable fact is that the above lower bound holds e.g. for 
compactly supported data (that is, even without assuming that \a{x)\ > c(l -|- |a;|)~'^~''^). In 
particular, this allows us to improve the previously known results (see e.g., [3], [19]) on the 
instantaneous spatial spreading property of highly localized flows. 

For generic global strong solutions. Theorem 1.7 implies various lower bounds. More 
precisely, starting from a fast decaying initial datum, we get, for all < a < d -|- 1, and 
large t : 

||n(t)||icj3 >ct-('^+^-"V2. (12) 

This result is a converse to Miyakawa's property (6). 

In the same spirit. Theorem 1.7 can be applied to estimate the decay of the moments of 
the solutions : for all 1 < p < oo and a > such that 

a + -<d+l, (13) 
P 
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we obtain, for large t, 

\\{l + \x\ru{t)\\L,>ct--^^''+'-''-i\ (14) 

This lower bound seemed to be known only in a few particular cases (namely, p = 2 and 
< a < 2, see [23], [12], [4], or 1 < p < oo and a = 0, see [8]). The corresponding upper 
bounds to (14), starting with the work of M. E. Schonbek, have been studied by many authors. 
See [18] for a quite general result. 

In some sense, the restriction (13) on the parameters could be removed, since (11) implies 
that for generic solutions, one has 

+ \x\ru{t)\\ L, = oo 

whenever a + ^ > d + 1 (see also [3]). 

Our results can also be applied to the study of the anisotropic decay of the velocity field. 
In the whole space, we show that not too stringent anisotropic assumptions on the decay of 
the data will be conserved by the flow. We also show that, if the initial data is well localized 
in M.^, then the flow decays faster than (1 + |x|)~^'^"'"''^^ as soon as one component does. This 
prevents localized flows in M*^ from having a really anisotropic decay. 

The situation can be different in other unbounded domains. For example, we will briefly 
discuss the case of the half plane > 0, with Neumann boundary conditions, and show that, 
in this case, generic flows have a genuinely anisotropic decay. 



The asymptotic separation of variables method 

The proof of (2) relies on a new, simple method that is a sort of "asymptotic separation of 
variables". We can summarize it as follows: one starts writing the Navier-Stokes equation in 
the usual integral form 

u{t) = e*^a - / e(*-^)^ Pdiv(u ® u){s) ds, (15) 
Jo 

where P is the Leray-Hopf projector onto the divergence-free vector fields : 

P/ = /- VA-Hdiv/). 
Then we use a classical decomposition of the nonlinear term (see e.g., [23], [3]) 

{ui^u){x,t) = { / {ui^ u){y,t)dy]g{x) + v{x,t), 

where q denotes the standard gaussian function, and v is defined through this formula. 

S,nce / = 0. the^nction .-Pdiv. behaves at i„iin.t, better than the ptevi- 

ous non-linearity e*^P div(;/ 5§ u) : its contribution can be included in the remainder terms. 
The next step consists in observing that the kernel of e*^P div behaves, as \x\ — >■ oo, as a 
time independent homogeneous tensor H{x). Then we show that by applying e(*~*)^P div 
to a matrix of the form S{s)g, where the coefficients of S{s) depend only on time, we get 
H{x) ■ £{s), plus some lower order terms. A time integration then yields a principal part 
for the velocity field of the form H{x) • K{t), as \x\ — oo. An explicit computation of this 
product provides the expression for Ii(x,t) in (2). 

We point out that the above strategy is not specific to the Navier-Stokes equations, but can 
be adapted to obtain the spatial asymptotics for more general models. What one essentially 
needs for its application are sufficiently explicit expressions (or sharp estimates) for the kernels 
of the operators involved. 
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Structure of the article 

Our main results are Theorem 1.2, its companion Theorem 1.7, and Theorem 3.1. Corollary 3.6 
also has some interest, since it extends a few results in the existing literature and its proof is 
very short. This paper is organized as follows. The proof of Theorem 1.2 and Theorem 1.7 
is contained in Sections 2.4-2.5, after we have prepared some preliminary estimates. In 
Section 2.6 we establish Proposition 1.6, in a slightly more complete form. The remaining 
part of the paper is devoted to applications : in Section 3.1 we give a precise statement and 
a proof of (11). Section 3.2 contains the proof of (12) and (14). The last sections deal with 
the anisotropic decay of solutions. 

2 Proof of the main results 

Let us now focus on the proof of the above results. 

We shall use the following notations for the kernel of the convolution operator e*^P div : 

According to (15), the j^^ component of (NS) can therefore be written as 

d d „t 

Ujit) = e^^CLj - / Pj;h,kit - s) * {uhUk){s) ds. (16) 

This kernel is related to the standard gaussian function g{x) = (47r)-<^/2e-l=^l'/'' in the 
following way. One has : 

with 

■j^ POO 

= 2 i(9h9t)Sj,k + {dk9t)hh] , -^S,fc(^' *) = djdhdkgs{x) ds 

Note that fJ.^^ k ~ ^j^k h ^j'^h k ~ ^j'^k h accordingly to the fact that only the symmetrical 
kernel has a physical meaning; gt{x) = t~'^/'^g{x/\/t) is the fundamental solution of the heat 
equation. 

2.1 Some elementary computations on F 

We shall need time-independent asymptotics of F, valid in the region where \x\'^ S> t. 

Lemma 2.1 There exist two positive constants C and c that depend only on d, and a family 
of smooth functions ^ j-^h,k satisfying 

\^j;h,k{^)\ + \^^r,h,k{^)\ < Ce-\-\' (17a) 

such that : 

PrAk{x,t) = 7. (^^^ - (^ + 2)|pij) + \-\-''^''^rAk{-^) 
with 7d = TT-'^/^ r(^^) and (Tj^h,k{x) = Sj^hXk + Sj,kXh + Sh,kXj. 
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Remark 2.2 Note that x = is not a singular value of F ; indeed, Fj.h,k is a C°° function on 
M'^x]0;+oo[ and one may immediately check on the Fourier transform that Fj-ii^k{0,t) = 0. 
Actually, for \x\'^ < t, the following computations also imply that : 

\Fr,h,k{x)\ < C -^-^^ + O (^^(^^ j ■ 

Remark 2.3 In Theorem 1.7, one has : 



Ej{y,t) = I I {uhUk){z,t - s)'^j.h,ki \ y] dzds 

with the ^ j-h,k given by this lemma. 

Proof. For all indices j, h, k (distinct or not) in {1, ... , d}, one has : 



F' 



:i) +\ _ ^j,kXh + SjjiXk i3:|2/4t 



thus F^^,ix,t) = |x|-(<^+i) ^(^^ [x/Vt), With 

Let us introduce (Tj^h,k{x) = Sj^hXk + Sj,kXh + 5h,kXj- One also has : 
^(2) u t) - r ( ^^'■^'"^'^^ - ^1^] a (x)ds 

The change of variable A = |x|/\/4s gives gg 

-dyi^ A ^ g^j^^ therefore 
The following formula provides information when A= |x|/\/t ^ 1 : 



f 

Jo 



This leads to 



^d/2p(2) .^ _ '^jAkjx) r.(d+'^\ 'iXjXhXk + I ..,(2) (x\ 

^ ^ r,h,k\^^^> - |^|d+2 ^ \^ 2 y \x\d'+^ \ 2 ^"^1 ' 
with 

' kl J\x\/2 \xr J\x\/2 

Conclusion (17b) follows immediately from the well known formula r{z + 1) = zT{z). The 
bounds on '^j;h,k and its derivatives are also obvious. 

□ 

The second valuable property of F is that the convolution with the standard gaussian 
function is equivalent to a shift in time. 
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Lemma 2.4 For all t > and x G M'', one has : 

{Fj-hM-^ t) * 9) (x) = Fj,nA^, t + 1). (18) 
Proof. Since {gt)t>o is a convolution semi-group, i.e. gt * g = gt+i, one has : 
FfU-^i) *9=\ [{dh9t)5i,k + {dk9t)5i,u] * 9 = f'^U; t + 1) 



/•oo roo 

and Ffl^{-,t) *g= djdhdkgs+iix) ds = djdhdkgs{x) ds = i^Sik(-,t + !)• 
' Jt Jt+i ' 



□ 



Let us finally recall a classical estimate of the norm of the kernel. 
Lemma 2.5 There exists a constant C > such that 

Vt>0, <CrV2. (19) 

Proof. This follows from (17b). 

□ 

2.2 Decomposition of the non-linear term 

Theorems 1.2 and 1.7 rely on a suitable decomposition of the non-linear term. A similar 
decomposition has been previously used by M. Schonbek [23] to prove lower bounds on the 
large-time decay of the L^-norm of the flow. This part of the computations is common to 
both proofs. 

Let us first explain the decomposition on a gaussian non-linearity. If ^ denotes the standard 
gaussian function and gt the fundamental solution of the heat equation, one sets : 

9hx) = (^j^^gt{yfdy^g + /\{x,t). 

The remainder A(x,t) has a mean value of zero : 

A{x,t)dx = 0. 



/ 



For fixed x G M*^, this approximation scheme behaves badly if i — > or i — > -|-oo ; indeed, a 
simple computation leads to 



A(x,t) = |l-Q)'^'eM^lV^*|5?(x). 

But this computation also shows that such approximation scheme is satisfactory at least when 

i~2 or IxP ~ 2o? — ^In^, 
' ' t-2 2 ' 

i.e. when A(x,t) is close to zero. We now perform a similar decomposition for the non- 
linearity in (NS). 
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Now let a be a divergence-free vector field and u G C^([0, r];L^) be the solution of 
(NS) given by Theorem 1.1, starting from these initial data. Recall that ■& = min{i?; d + 1}. 
Since > ^, one has C and hence the energy matrix 

is well defined. Consistently with the preceding approximation scheme, let us define Vh,k by 

{uhUk){x,t) = Sh,k{t)9{x) +Vh,k{x,t). (20) 
Thanks to Lemma 2.4, the integral equation (16) is hence equivalent to : 

Uj{t) = e^^aj - y] / £h,k{s) Fj.h,k{^ + 1 - s) - V / Vh^k{s) * Fj;h,k{t - s) ds. 
h,k h,k 

The time-independent asymptotic (17b) of the kernel Fj.h,k now leads to 

\x 

where Pj is given by 



Uj{x,t) = e^^ajix) + -^j^ +Rjix,t), (21a) 



Pj{x,t) = 7dX] {id + 2)xjXhXk - \xfaj^h,ki^)) ^h.kit), (21b) 
h,k 

with aj^h,k{x) = 6j^hXk + Sj,kXh + h,kXj. The remainder 

R^{x,t) = 



h.k 



is given by : 

= kr^'+'^ f^£h,k{t - ^) *^;'^.^(^=y) ds, (21c) 

^S,fe(^'*)= fvh,k{s)*Fj.f,,k{t-s)ds. (21d) 
Jo 

The functions '^j;h,k are given by (17b). 

Remark 2.6 The above remainder is not small when \x\ < As the solution u{x,t) is 
smooth at least for small t > 0, the homogeneous polynomial and the remainder have to 
behave in exactly anti-symmetrical ways when |x| — > 0. The same compensation also occurs 
for a.e. x eW'' when t — 0. 

The polynomial profile P{x,t) = \x\~'^~^P{x,t) has no vorticity, i.e. the matrix 

rotP=(a,p._a.p^).. 

is identically zero. This means that the polynomial profile is a gradient vector field. In fact, 
one may check immediately that : 

^ = vn with n(x,*) = 7i^-Eraa-^''.M*))- P2) 
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2.3 General bounds of the remainder terms Rj-h,k 

Let us now compute some upper bounds of the remainder terms. This second part of the 
proof is also shared by Theorem 1.2 and 1.7. 



Bound of Rj^^k{x,t). The bound (17a) gives 



^j h k — , < C exp , ^ , , , , 



I 1 2 

C \x\ 



hence 



Bound oiR):^l^^^{x,t). Since Vh^kip^i s) dx — 0, the second remainder can also be written : 
^?hki^'^) = f f Vh,k{y,s) {Fj-h,kix -y,t-s)- Fj.h,k{x,t- s)) dsdy. 

■" ' Jo JRd 



The Taylor formula gives : 



l-^?iik(^'*)l ^ / / \y\\'Vh,k{y'^)\ sup |V-Fj;/i,jt(x + 2;,t-s)|ds(iy 

•''0 J\y\<\x\/2 \z\<\x\/2 

+ f if \vh,k{y^s)\dy]\Fj.h,k{^,t- s)\ds (24) 

Jo \J\y\>\x\/2 J 

+ [ [ \vh,k{y,s)\\Fj.h,k{x-y,t- s)\dsdy. 

Jo J\y\>\x\/2 

Thanks to (17a)-(17b), one has \V Fj.h,ki.x,t)\ < C\x\-^'^+^'^ uniformly for t > 0. Applying (19) 
as well, we get : 



\Rflki^M<c( f f \y\\vh,k{y,s)\dyd^\x\ 

■'' ' \JoJ\y\<\x\/2 J 



-(d+2) 

\y\ \'^ii,K\yi ■^n j i-^i 

'\y\<\x\/2 

+ ( [ [ \vh,k{y^s)\dyds]\x\-^''^'^ (25) 

\7o^|j/|>|a;|/2 / 

+ / it-s)~^^'^ sup \vh,k{y,s)\ds. 
Jo \y\>\x\/2 

To conclude the proofs of Theorem 1.2 and 1.7, we shall now use the assumptions on u to 
estimate (23) and (25). 



2.4 Local-in-time solutions. Proof of Theorem 1.2 

The goal of this section is to get upper bounds of the above remainders that provide valuable 
information for short time. In particular, in view of the proof of the lower bounds (11), it 
is of interest to have information on the behavior as t — ^ of the last term appearing in the 
right-hand side of (2). 

The remainder satisfies 

^ ds. (26) 
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Indeed, if < < one has 



|2 

/ L. I 

exp 



and if ?? > , one also has 



.1/2 



In both cases, our estimates can blow up as i ^ oo, but not faster than (1 + 1)^/'^ . 

(2) 

To deal with the remainder Rj.j^y., one may notice that the definition of v implies 

\vh,k{y^^)\ < + Hs)\\l9iy) < c(i + lyl)-^^'' 11^^(^)111^= 

with ?? = min{i? ; d+ 1}. Therefore, since 2'& > d+1 : 

\4lk{x,t)\<C\x\-^-' f\\u{s)\\l^ ds 

J 

+ + f\\u{s)\\' {t-sry^ds, 

Jo ^ 



and hence for \x\ > 1 



< -p . (27) 



This ends the proof of (2). To obtain an asymptotic profile for the pressure, we need the 
following simple result on the localization of the derivatives. 

Proposition 2.7 Given u G L°°([0, T]; a solution of the Navier-Stokes system with 
Cauchy datuna a = u{0) and < i} < d + 1. If, for some index i, one has dia £ L^, then : 

diueL'^{[0,T];Lf). (28) 

If, moreover, dia G and didja G L'^ holds for all j G {1, . . . , d}, then 

didjueL^{[0,T];Lf) (29) 

and 

t^/^dtueL^i[0,T];Lf). (30) 

Proof. Let us first deal with the first order spatial derivatives. Taking the i-th derivative 
in (16) leads to the aflfine fixed point problem : 

diU = G{diu) 

with e = {Qi,...,Qd) and 

= e*^{diaj) - 2 V / -P,;/i,fe(i - s) * {uhWk){s) ds. 
h,k 
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Proposition 3 of [24] implies that 6 is a continuous operator on X = L°°([0, Tq]; L^), < 
To < T and that 

\\e{w - w')\\ < CoT^ sup ||M(i)||ioc ||u;-i(;'|L. 

te[o,T] ^ 

One may therefore choose Tq > such that is a contraction of the Banach space X. Its 
only fixed point w = diu belongs therefore to this function space. The same argument also 
holds on [To, 2To],. . . and leads finally to diU G L°°([0, T]; L^). 

Conclusion (29) also follows from the contraction mapping theorem in a similar way. 

For the time derivative, the starting point is again an identity that directly follows 
from (16), namely, 

dfU = Q{dtu), 

with 

@j{w) = Aaj - ^ Fj.h^k{t) * (ahak) - 2 ^ / Fj.h,k{s) * {uhWk){t - s) ds. 
h,k h,k •'^ 

The Banach space we deal with is 

y = {«;;ti/2|mt)|| GL-([0,ro])}. 
Proposition 3 of [24] implies now that 

-Aoj + ^i^i;/i,fe(t) * (a/iflfc) e Y 
h,k 

and 

1 

||0(i(; — 'ii;')||„ < ttCo Tq^ sup ||ti(t)|| roo — . 

te[o,T] 



Vi > 0, [ —^^^= = TT. 

Jo 



Here, we have used the fact that : 

ds 

^/s{t-s) 

The conclusion now follows in the same lines as above. □ 
We can now establish (4) : The pressure is defined up to an arbitrary function of t by : 

—Vp = {dt — A)u + div(M (Si u). 

Let us now replace u by the its profile given by (21a), that \s u = e^^a + VII + R. One gets : 

-Vp = V{dt - A)n + {dt - A)R + div(n (g) u). 

This yields : 

-J^d- - • ^h,k{t) + t), (31) 

h,k ' ' J 

where the remainder term q{x,t) satisfies 

-Vq = -VAU + {dt- A)R + div{u^u). (32) 
Let us show that, for all i > 0, we have Vq = Ot(|a;|- ™°{2i9,(i+2})^ 
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- We obviously have VAH = min{2i9,d+2}^^ since the left hand side is a homoge- 
neous function of degree —d — 3 which is smooth for x ^ 0. 

- The term div(M u) = {u • V)u belongs uniformly to because of (28). 

- The remainder is the sum of two terms: one checks immediately that {dt — A)R^^^ is 
exponentially decaying as \x\ — >■ oo. The second term is 

{dt - ^)Rfl^k^x, t) = Vh,m * Fj.h^kit) + J^i^t- ^)vh,k{t - s) * Fj,h,kis) ds 

where Vh,k{x,t) = uuUk - £h,k{t) g{x) and t^/'^dtVh,k belongs to L°°([0, T]; L^). 

We now use again that j Vh,k{^) dx = : if we apply Lemma 2.5, the computation (25) 

shows that the first term is bounded in i~^^^-^°°([0, T]; L^jj^|2^.^_l_2j.). On the other 
hand, dfVfi^k and Avh,k also have a vanishing integral. Therefore, using the estimates 
on the space-time derivatives provided by Proposition 2.7 shows that the second term 
belongs to i^°°([0,r];L~„^2^.,^2})- 

Hence we get Vg G Ot{\x\~™^^^^'^''^'^'^^) . Our last step is the following elementary estimate : 

Lemma 2.8 Let a > 1 and f G C^{W^) such that V/ G L'^ . Then there is a constant c such 
that f - c e L'^_i. 

POO 

Proof. For any u G M*^, |a;| = 1, let = limj._»oo f{f'^) = /(O) + / V/(sa;) • ujds. If ui is 

Jo 

another point of the unit sphere then for all r > we have 

poo f'OO 

\L-ico\< \Vfisio)\ds + Cr sup \Vf{x)\+ \Vfisu)\ds. 

Jr \x\>r Jr 

Letting r — > oo we get that c = is independent of to. But 

\-a+l 



/•oo 

Ifiru) -c\< |V/(sa;)| ds < C(l + r)' 

Jr 



and the conclusion follows. □ 

The standard properties of strong solutions imply that q{x,t) is smooth for x ^ 0. Ap- 
plying this lemma (for fixed t, < t < T), with f{x) = x{x)q{x,t), where X is a smooth 
function such that x(a^) = for |x| < r and x = 1 for > r' for some < r < r' implies 
that, g(a;,t) = c-FOt(|x|-™°{2'^-^''^+^>). This completes the proof of Theorem 1.2. □ 

For later use, let us note explicitly that if R{x,t) = Ot(|x|~™'''{2'^''^+^>) denotes the last 
term in the right hand side of (2), then we have proved that, for all < t < T, 

C(Vi + t) \\a\\Lr 
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2.5 Global-in-time solutions. Proof of Theorem 1.7 

Let us now focus on long time asymptotics. Let u be a global solution satisfying (6). Going 
back to (21a), we see that the profile (7) holds with 



{y,t) = - f I i'^hUk){y, t - s) ■^j.h,k { \ y] ^^^^ 

and 

h,k 

The bound (17a) immediately implies (8). To prove (9), we start by observing that 
assumption (6) implies : 

\Ht)\\l,<c,{i + t)-i^-^-') 

for any e > and therefore, letting v = (vh^k), 

\viy, s)\ < \uiy, s)|2 + \\uis)\\l g{v) < C {I + \y\r^'^ (1 + s)-^^'-) 
for d/2 < a < i?. 

A consequence of (17b) is that, for all < /? < d + 1 and < 7 < 1 

|F,;,,fc(x,^)|<c|xr'^^^('^+l-w^ 

|VF,;,,,(a;,t)|<|x|-('^+^+^)i-M/^ 
Therefore, coming back to (24) we get 

\n{x,t)\ <c( f ( (l + |?/|)i-2°(l + s)-''+"(t-s)-5+l ci?/dsVr^'^^^+^^ 

\jQj\y\<\x\/2 J 

+ c( [ [ il + \y\)-^''il + s)-^+'^\x\-^{t-s)-^dyds) 

\JoJ\y\>\x\/2 / 
+ C [\t - s)-^/^(l + k|)-2°^(l + 5)-^^+" ds, 

Jo 

where, in the last integral, we have also used Lemma 2.5. 

Let us call Ii, I2 and 1^ the three terms of the right-hand side. To estimate Ii, we fix a 
small e > and choose a = + e. Then we write Ii = Ii^i + Ii^2, where these two terms 
are obtained by splitting the integral into jj^^ and J^^^- Then we have, for alH > 1, 



(34) 



and 
Thus, 



h<C,,e\x\ <!^^_,+,^ if^<^<^, ^^^^ 



To estimate I2, we choose again a = + e and P = d + j. Then the same argument as 
before shows that I2 can be bounded as in (35). To estimate I3, we take a = '^"'"^"'"^ ■ This 
choice shows that Is is also bounded by the function on the right hand side of (35). Summing 
all these bounds completes the proof of Theorem 1.7. 

□ 
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2.6 Criterion for the vanishing of VH - Proof of Proposition 1.6 

We now prove Proposition 1.6, which we restate in a more complete form. 

Proposition 2.9 For any real matrix K = (Kfi^k), l^t us define a family of homogeneous 
polynomials by 

Qji^) = XI -{d + 2)xjXhXk) Kh,k- (36) 

h,k 

The following assertions are equivalent : 

1. The matrix is proportional to the identity matrix, i.e. 

\/h,k £ {l,...,d}, Kh,k = aSh,k (37) 

with a = ^Ti K. 

2. Qj = for all indices j £ {1, . . . ,d}. 

3. There exists an index j € {1, . . . ,d} such that Qj = 0. 

4. There exists an index j E {1, . . . ,d} such that djQj = 0. 

Putting the terms factor in (36), one gets the following expression for the j™ 

component of Q : 

d 

Qj{x) = Xj J] {Tr K - dKe,e + 2{Kjj - Ke^e)} xj 

+2\xfk{Ej,x) -{d + 2)xjk{x, x) 

where Ej = (0, . . . , 0, 1, 0, . . . , 0) denotes the canonical basis of and K is the bilinear form 
defined by the non-diagonal coefficients of K : 

K{u,v) =^Kh,kUhVk- 
h^k 

Relations (37) express the fact that the matrix K = {Kh^k)i<h,k<d is a scalar multiple of 
the identity matrix. In such a case, one can immediately check on the previous expression 
that Qj{x) = 0. 

Let us prove conversly that djQj = implies = aid. One has : 

d 

djQj{x) = X + 2<^J-,')(T> K - dKe,e) + 2{Kjj - Ke,e)} xj 

i=\ 

- 2dxjk{Ej,x) - {d + 2)K{x,x). 
The fact that djQj{Ki) = for all i implies 

G {1, . . . , d}, (1 + 2dj,i){Tr K - dKe,e) + 2{Kjj - Ke,e) = 
and hence Ki^i = ^Tr K {i = 1, . . . ,d), i.e. all the diagonal entries of K are equal. Therefore : 

djQj{x) = -2dxjK(Kj, x)-{d + 2)K{x, x) 
and this expression should vanish identically. A new derivation with respect to Xj gives 

d'^^Qj = -A{d - l)k{Ej,x) = 0, 

i.e. K(Ej,x) = as d > 2, and hence K{x,x) = 0. This proves that the matrix K is a scalar 
multiple of the identity matrix. 

□ 
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3 Applications 



or such that 



Let us now explore a few consequences of the above results. 
3.1 Instantaneous spreading property 

It is a consequence of the result of [3] that, if the components of the initial data have no 
special symmetries, then the corresponding solution u{x, t) of (NS) satisfies 

hminfi? / \u{x,t)\dx > , (38) 

Jr<\x\<2R 

for t > belonging at least to a sequence of points tk converging to zero as — > oo. In 
particular for those t, one has : 

/ NM.,t)l<fa = +=c and / rf«Wx.*)l^d. = +oo. 

(See also, e.g., [19], Theorem 25.2). The precise assumption guaranteeing (38) is the non- 
orthogonality of the components with respect to the L^-inner product, i.e. one can find j ^ k 
in {1, ... , d} such that 

/ aj{x) ak{x) dx ^ (39a) 

/ a]{x)dx^ / al{x)dx. (39b) 

Even if (38) already explains that the limitation '& < d + 1 in Theorem 1.1 is optimal for 
generic flows, such a condition does not provide much information on the pointwise decay of 
u, as \x\ — ^ oo. 

Theorem 1.2 and Proposition 1.6 not only provide such information, but also allow us to 
give a simpler proof of these facts. The instantaneous spreading property is fully described 
by the following result. 

Theorem 3.1 For > d + 1, let a G be a divergence-free vector Geld. Let u be the 
corresponding solution of (NS) in ([0, T]; L^^) . For <t<T, we set 

«:t = max{l,rV2,t-V('?-<i-i)}. 

1. There is a constant c > such that, for < t < T and \x\ > CKf : 

\u{x,t)\<ct\x\-^'^+^\ (40) 

2. Conversely, if (39) holds for a couple of indices {j, k), then there exists to G (0, T] and a 
constant d > such that for all < t <to and all x in a conic neighborhood of the Xj 
or Xk axis, with \x\ > CKt : 

\uj{x,t)\>c't\x\-'^'^+^l (41) 

3. Actually, if (39) holds, the lower bound (41) holds in almost all directions : the set 

(42) 



a e S'^-^ ; liminf (t'^ \x\'^+^ \uj{x,t)\J = 

+0O 
x&La 

is a closed subset of the sphere , of measure zero. 
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Remark 3.2 The non-orthogonality assumption (39) cannot be removed. Indeed, [3] contains 
an example of symmetric and highly localized flow for which (38) and (41) break down. For a 
further understanding of how symmetries usually lead to highly localized flows, one may also 
refer to [16] or [2]. 

Tlie previous results show that, for generic flows, each component of the velocity fleld 
instantaneously spreads out in any direction. In particular, one has the following corollary. 

Corollary 3.3 Under the same assumptions as in Theorem 3.1 there exist a time to > and 
two positive constants cq, C such that for all t G {0,to] : 



Cot 

\Uj[x,i) \ ax ^ 

IR<\x\<2R 

for all R > ct-i/min{2;,?-(i-i} particular, for j,k = l,...,d 



L 



Uj{x,t)\dx>-^ (43) 
R 



I 



Xk f P\uj{x,t)\P dx = +00, (44) 



as soon as 1 < p < +oo and '& + ^ > d + 1. 

Remark 3.4 The spreading property (44) should be compared with F. Vigneron's result [24], 
according to which 

/ {l + \x\fP\u{x,t)\Pdx <Ct [ {l + \x\fP\a{x)\Pdx 



as long as the solution exists in /^(R ), provided that p > d and + ^ < d+1. 

Proof. First, let us establish the upper bound (40). Theorem 1.2, together with (33), imme- 
diately implies, for |rc| > 1 : 

\u{x,t) - e*^a{x)\ < Ci|a;|-('^+i) < Ct\x\-^'^+^\ 

for all xeW^ such that |x| > i^^^. Moreover, if |x| > (2i-i)V('^-'i-i) ^ one has : 

|e*^a(x)| < C|a;|-^ < (45) 

This proves (40). Let us now focus on the lower bound (41) of \uj{x,t)\. 

Let j ^ k such that a = J^a{ajak) dx ^ 0. Then, for some to > 0, possibly depending 
on a, and all < t < to, we have 

\Kim > 

Let e > and let = {x; \xr\ < e\xk\ (r / k)} be a conical neighborhood of the x^-axis. 
Recalling (21b) we get, for e small and R large enough : 

VxGFfe, \x\>R =^ \Pj{x,t)\>^-^\xkf. 
Using (21a), (33) and the flrst of (45) now leads, for large enough |x| and x G F^, to 

\u{x,t)\ > \uj{x,t)\ > > \^\x\-id+i) 

and (41) follows in this case. 
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In the second case we choose j ^ k such that j^^i <i] dx 7^ Jj^^ dx. Then we have 

d dx ^ y / a^dx 

(otherwise, d j a^dx ^ Ylm=i /r^ dx, and we should exchange j with k). Let us set 

(5 = \^ / a^dx — d a] dx. 



Arguing as before, we see that there exists a conic neighborhood Tj of the Xj axis, such that 
for all X G Tj and |x| large enough we have 

\u{x,t)\ > K-(x,t)| > M!|^.|-(<i+i) > M!|^|-(<i+i). 

Then (41) follows in this second case as well. 

Let us now prove the last statement of Theorem 3.1. The map s {uhUk){x, s) dx is 

continuous. Therefore, (21b) implies that 

Vx G R"^, lim -P(x, t) = Vix) 
t^o t 

where V = {Vi , ■ ■ ■ Vd) is given by 

Vjix) = Td ( / ahttk ] {{d + 2)xjXhXk - \x\'^aj^h,k{x)) , 

which is a homogeneous polynomial of degree exactly three. According to Proposition 1.6, 
the assumption (39) means that Vj ^ for all j = 1, . . . , d. 

The convergence of P{x,t) to V{x) is uniform when x belongs to the (compact) unit 
sphere. Let us define the following dense open subsets of S"'"^ : 

nj = {ue S'^-^: Vj{uj) ^ 0}. 

Given w G Q^, let us define T^, > as the supremum oit<T' such that 

Also let = \\Vj{uj)\. 

Prom (21a), (33) and the obvious estimate |e*^a(x)| < C|a;|~'', we get, for u = x/\x\ G Q,j : 

\uj{,x,t)\ > 2c^t |xr(<^+i) - CtV2|^|-{<i+2) _ c\x\-'^ > c^t Ixp^'^+i) 

for all < i < and > C^t~'^* . The complement of VLj is an algebraic surface and 
therefore the set S = has measure zero in S*^"-*^. 

Finally, (43) and the corollary follow immediately from (41) and the fact that a function 
bounded from below, at infinity, by Ixl""^"^ does not belong to any weighted Lebesgue space 
i7(M'^, (1 + |x|)'^*' dx) when l<p<ooandi? + ^>d+l. 

□ 
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Remark 3.5 The set of exceptional directions where, at some time t > 0, our lower bound (41) 
breaks down corresponds to the roots of the polynomial P{-,t) on the sphere In dimen- 

sion two we have a more precise description of such a set : it is made of at most six directions . 
This is due to the fact that a homogeneous polynomial in of degree exactly three has at 
most six zeros on S^, as is easily checked by passing to polar coordinates. In principle, this 
conclusion is valid only during a short time interval [0,to]- Indeed we cannot exclude that 
after some time there are flows featuring some kind of creation of symmetry: it may happen 
that, for some time ti > to, f ui{x,ti)'^ dx = f U2{x,ti)'^ dx and [{uiU2){x,ti) dx = 0, in a 
such way that the velocity field first instantaneously spreads out but then recovers a good 
localization at time ti. Subsequently, the flow would remain localized, or spread out again, 
depending on whether the components ui and U2 remain orthogonal or not, after ti. However, 
no example of such a somewhat pathological flow is known so far. 



3.2 Lower bounds of solutions in weighted spaces 

Let us establish a few consequences of Theorem 1.7. Throughout this section we suppose 
t>l. For well localized data, e.g., when a G Lf with i? > the limit 

Vnoo(x) = lim VU{x,t) 

t— »oo 

is well defined. In this case, a consequence of (7) is that, for some P > 0, 

< C|x|-'^-ie-l"l'/(*+i) + C\x\-^-' t-^. (46) 
This leads us to introduce, for all ^ > 0, the region 

VA{t) = {xeR'^; \x\'^ > A{t + 1)}. 

Since, generically, VIIoo ^ 0, several lower bounds for the large time behavior of u can be 
obtained as an easy consequence of (46). For example, if we introduce the weighted norm 

ii/iils = (/ \f{x)ni + wdx 

then, taking A> large enough, for all 1 < p < oo and a >0, such that 

a + -<d+l, (47) 
P 

we get, for alH > large enough. 




> ^ / \VU^{x)\P{l + \x\r''dx (48) 

> c {At)-'^^''^'-''-i\ 

If the datum is highly oscillating (for example, if the Fourier transform of a satisfies some 
suitable vanishing condition at the origin) then the norms of e*^a decay faster as t ^ oo 
than the right hand side of (48). Then, (48) will be in fact a lower bound for ||n(t)||% in 
this case. A similar conclusion remains true if we drop this assumption on the oscillations 
and start with a datum that is simply well localized. Indeed, we have the following : 



u{x,t) - e*^a(x) - Vnoo(x) 
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Corollary 3.6 Let u he as in Theorem 1.7, starting from a G L"^, with i? > d + 1. We 
also assume that VHoo ^ 0. Then there exist to > and a constant c > such that, for all 
1 < p < oo and a > 0, satisfying (47), we have, for all t>to: 

lhWllLS>cr^('^+^-"-?). (49) 

Moreover, for allO < a < d+ 1 and all t > to : 

lk(t)||i^ >ci-5(<^+i-"). (50) 

The lower bound (49) was already known for p = 2 and < a < 2 (see, e.g., [23], [4]), 
or 1 < p < oo and a = (see [8]). The decay profile (6), under the assumption of Corollary 3.6, 

immediately implies the (slightly weaker) upper bound ||«(t)||iP < c^r^^'^+^-^^-i;-^) for 
all e > 0. In fact, the "sharp" upper bound {i.e. the bound with e = 0) has been obtained, at 
least for p > 2 and with some additional restrictions on a, by many authors (see [18] and the 
references therein). 

Proof. By our assumptions, a G L^{W^) and diva = 0. Thus, J a{y) dy = 0. A direct 
computation (using the same method as in the proof of (24)) then yields 

|e*^a(x)| < C(l + |a;|)-^ (1 + t)('?-'i-i)/2. 
It then follows that, for t > 1, 

/ |e*^a(x)r(l + \x\rdx < 

Here the exponent of A is strictly smaller than that of (48). If A is large enough, then a 
comparison between this inequality and (48) gives (49). The proof of (50) is essentially the 
same. □ 



3.3 Flows with anisotropic decay in the whole space 

This short section contains a positive and a negative result about flows in with anisotropic 
decay at infinity. 

Theorem 1.2 implies that (NS) flows may inherit the anisotropic decay properties of the 
initial data, as long as these properties do not violate the instantaneous spreading limit given 
by Theorem 3.1. 

Proposition 3.7 Let a be a bounded divergence-free vector field and u the corresponding 
solution of (NS) in ([0,T);L°°). Let us also assume that there exists a function m such 
that 

\e'^a{x)\ < Ct{l + \x\)-^m{x)-^ (51) 

with ^<-&<d+landl< m{x) < C(l + \x\f+'^-^. Then, for aU T' < T there exists a 
constant Ct' (this also might depend on the data) such that : 

\u{x, t)\ < Ct'{1 + {xl)-"^ m{x)-^ (52) 

for all te [0;T']. 
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Proof. This is an obvious consequence of (2) and (33). 

□ 

Let us give some examples of anisotropic weights satisfying (51). A Peetre-type weight is 
a measurable function m : R'' — [1; +oo) such that 

3Co > 0, Vx, y e R'^, m{x + y) < Co m{x)m{y). (53) 

Common examples are (for > 0) : 

mi{x) = 1 + Ixil""^ + . . . + \xd\'^^ and m2(x) = e^^l^L 

The class of Peetre-type weights is stable by finite sums and products, translations and or- 
thogonal transforms. 

Lemma 3.8 Let m be a Peetre-type weight such that m{x) < C exp(c|x|) and T > 0. Then, 
there is a constant Ct > such that 

||m(e*^a)||^^ <CT||ma||ioc. (54) 

Proof. It is an elementary computation : 

m{x) |e*^a(x)| < Co [{mgt) * (m|a|)] (x) 

< Co (47r)-'^/2 II 

The conclusion follows from the bound m(^\/iy) < Cexp(cr|y|). 

□ 

As a converse to the previous result, the following property implies that highly localized 
flows cannot decay at inflnity in a really anisotropic way. 

Proposition 3.9 Let a £ -L^^+e ^ divergence-free vector held with < e < 1, and u the 
corresponding solution of (NS) in ([0, T); L^j^) . For some t > 0, let us assume that there 
exist an index _7 G {1, . . . , d} and a subset S C of positive measure such that 

Vo-gE, lim \xf+^Uj{x,t) = 0. (55) 

|3;|^+oo 

Then, there exists a constant C > such that 

K(x,t)|<C(l + H)-(<^+i+^) (56) 

for all k = 1, . . . ,d. Moreover, if (55) holds for a finite time interval t G [Tq,Ti], then C may 
be chosen uniformly with respect to t. 

Proof. Our assumptions imply that the polynomial Pj{x,t) identically vanishes. Proposi- 
tion 1.6 then implies that all the other components of P(x,t) also vanish. Our statement is 
once again a consequence of (2). 

□ 
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3.4 Application to the decay in a half-space domain 

Our last application of Theorem 1.2 is the study of the decay of solutions of the Navier-Stokes 
equations in the half space 

= {{x',Xd) : x' e Xd > 0}. 

We set u' = (ui, . . . , tx^-i) and x' = (xi, . . . , Xd-i). Let {e~^^'}t>o be the semigroup generated 
by —A' = A, in the case of the Neumann boundary conditions : 

ddu'\Q^d_ = 0, Ud\Q^d_ = 0. (57) 
where dd = The integral formulation of the Navier-Stokes system in is 

u{t) = e-^^'a - f e-(*-")^'Pdiv(tx u){s) ds, (58) 
Jo 

with diva = 0. We refer to [9] for the construction of weak and strong solutions to (58). 

We have the following result: 

Proposition 3.10 Assume that a € LfiRf), with ■& > Tien there exist T > and a 

unique strong solution u E Cw ^[0, T); L^(R^)^ of (58). Such a solution satisfies 

u{x, t) = e*^a(x) + H{x, t) + Ot (\x\- '^^^'^'^ ' '^+2>) , (59) 

where H = (Hi, . . . , Hd) is homogeneous of degree — (d + 1) for all t G [0, T), and such that: 

\Hj{x,t)\ < C\x'\-\x\-^'^+^\ {l<j<d-l) (60) 
\Hdix,t)\ < C\xd\ ■ \x\-(^+^l (61) 

As an immediate consequence of Proposition 3.10, we obtain the following anisotropic 
decay estimates (assuming that a is well localized). 

Ud{x,t) = C'(|x'|-('^+2)), u'{x,t) = C'(|x'|-('^+i)), when \x'\ +oo, Xd fixed 
Ud{x,t) = 0(|a;d|"('^+^)), u'{x,t) = 0(|xrf|-('^+2)), when Xd +oo, x' fixed. 

It is also worth noticing that Proposition 3.9 is not violated as the above decay holds only in 
a cylindrical region, and not in a conical one. 

Proof. This is immediate. Indeed, the study of (58) is reduced to that of (NS) in the following 
way. If u solves (58), then one can construct a solution of (NS) in the whole W^, setting 

Uj{xi, . . . ,Xd-l,-Xd,t) = Uj{xi, . . .,Xd-l,Xd,t) 

for j = 1,... ,d - 1 and Udixi,...,Xd-i,-Xd,t) = -Ud{xi, . . . ,Xd-i,Xd,t) (see [9]). Then 
under the assumptions of Proposition 3.10 we can apply (2) to u. But the integrals Kj^d(t) = 
jQj-gd{ujUd){x,s)dxds vanish, for j / d. Hence, from (21b) we see that \P{x,t)\ is bounded 
by a function H{x,t) satisfying (60). 

□ 
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Conclusions 



Theorem 1.2 provides a quite complete answer to the spatial decay problem of solutions of 
the free Navier-Stokes equations in the whole space, at least for well localized data. It would 
be interesting to know if some of the results of the present paper can be adapted to flows in 
other domains. 

For example, in the half-space case, the Neumann boundary condition considered in the 
previous section is not the most interesting one, since it destroys the boundary layer effects. 
The construction of the asymptotics as in Theorem 1.2, in the case of Dirichlet boundary con- 
ditions, would require a careful analysis of Ukai's formula (or its more recent reformulations) 
for the Stokes semigroup. 

In the case of stationary flows, asymptotic profiles have been given, e.g., by F. Haldi and 
P. Wittwer [15], [26]. Their results model the wake flow beyond an obstacle. However, they 
do not deal with the obstacle itself, but with a half-plane domain and a technical boundary 
condition dictated by experimental knowledge. 

For the non-stationary equation (NS) in M^\r2 with Dirichlet boundary conditions on 90, 

it seems reasonable to expect that anisotropic lower bound estimates for the decay of u should 
hold, when the net forces exerted by the fluid on the boundary, i. e. 



(where Tj^k[u,p\ = djUk + dkUj — Sj^kP and T\u,p\ — {Tj^k[u,p\)j,k is the stress tensor) do not 
vanish. This last condition, which is motivated by the results of Y. Kozono [17] and C. He, 
T. Miyakawa [16] on the L^-summability of solutions, would play, in the exterior domain case, 
a role equivalent to the non-vanishing criterion given by Proposition 1.6. 
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